CD 

ro 

ro 

co 

CQ 

O 

03 

. oa 



O 

< 




NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2902 



MATRIX METHODS FOR DETERMINING THE LONGITUDINAL -STABILITY 

DERIVATIVES OF AN AIRPLANE FROM TRANSIENT FLIGHT DATA 

By James J. Donegan 

Langley Aeronautical Laboratory 
Langley Field, Va. 



TECH LIBRARY KAFB, NM 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECH LIBRA"* 





n6bS*ttv» 


TECHNICAL NOTE 2902 


MATRIX METHODS FOR' DETERMINING THE LONGITUDINAL-STABILITY 
DERIVATIVES ’ OF AN AIRPLANE FROM TRANSIENT FLIGHT DATA 
By James J. Donegan 


SUMMARY 


Three matrix methods are presented for determining the longitudinal- 
stahility derivatives from transient flight data. One method, which 
requires four measurements in time-history form and utilizes the incre- 
mental tail load to separate the pitching-moment derivatives C m g and 

Cm^., permits the computation of all the longitudinal-stability deriva- 
tives. A second method requires three, measurements and one supplemental 

G m. 

assumption, namely = Constant. This method gives the most infor- 

ms 

mation for the least amount of work. The third method requires two 
measurements and two supplemental assumptions, namely = Constant 

and C m ^ = ^ Cpg (where Cmg and Cp^’ are the elevator-effectiveness 

derivatives, x^. is the tail length, and c" is the mean aerodynamic 

chord) . An Inspection of the results obtained, for the various methods 
shows the scatter which is typical of this type of analysis of flight 
data. 


INTRODUCTION 


The determination of the longitudinal-stability derivatives from 
flight data is a relatively difficult task because the wind-tunnel 
technique of permitting only one variable to change at a time, while 
constraining all the rest of the variables, cannot always be used. It 
is in the analysis of such flight-test data that matrix techniques 
employing the equations of motion seem to be particularly useful. 

Currently, much work is being carried out on the determination of 
stability derivatives directly from flight data but as yet this work is 
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still in the preliminary stages. The matrix methods for the • determina- 
tion of stability derivatives from transient flight data that are 
developed herein are an addition to this work. The previous work done 
on the determination of longitudinal-stability derivatives is extensive, 
and no attempt is made to summarize it since this summarization has beep 
adequately done in reference 1. 

In the present paper three methods are developed and presented for 
determining the longitudinal-stability derivatives from transient flight 
data. In these methods the expressions for some of the stability deriv- 
atives are in the form generally used in stability calculations. The 
first method requires the combination of four measurements in time- 
history form, two of which must he incremental elevator deflection and 
incremental tail load and the other two measurements can he chosen from 
a possible three, namely incremental load factor, pitching velocity, 
and angle of attack. The method demonstrates the use of the tail load 
to separate the pitching-moment derivatives and and deter- 

mine the - downwash derivative de/da. 

The second method, which is more restricted, requires a combination 
of three measurements (in time-history form), one of which must he incre- 
mental elevator deflection and the other two measurements can be chosen 
from a possible three, namely incremental load factor, pitching velocity, 
and angle of attack. This method also requires one supplementary assump- 
tion, namely = X C m ^, where X is a constant. 


The third method uses a combination of two measurements ( in time- 
history form), one of which must be incremental elevator deflection and 
the other one may be chosen from incremental load factor, pitching 
velocity, angle of attack, and so forth. The method also requires two 


supplementary assumptions, namely = XC^ and C m g = ^ (where 



and 



axe the elevator-effectiveness derivatives. 


X(- is the 


tail length, and cf is the m&an aerodynamic chord). By using a modi- 
fication of the third method, it is shown that considerable information 
can be obtained from a single time history. 


The methods are demonstrated by applying them to flight data 
obtained from tests of a medium jet bomber, and a comparison of the 
derivatives obtained by the various methods gives an indication of the 
accuracy which can be expected from data analysis by matrix techniques 
based on the longitudinal equations of motion. 
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b 

c 

C 1,C 2 

c L 

% 

% 

CLe 

Cm 

Cing 

Cmg 

Cn \i 

Cm t 

g 

I 

V 

L 


wing span, ft 
mean aerodynamic chord 
constants defined in appendix E 
lift coefficient, b/qS 

rate of change of airplane lift coefficient with angle of 
attack per radian; see appendix E 

rate of change of lift coefficient with elevator deflection 
per radian; see appendix E 

rate of change of lift coefficient with & per radian; see 
appendix E 

rate of change of lift coefficient with pitching velocity 
per .radian; see appendix E 

pitching-moment coefficient of airplane, M/qSc 

rate of change of pitching-moment coefficient with angle 
of attack per radian; see appendix E 

rate of change of pitching-moment coefficient with elevator 
deflection per radian; see appendix E 

rate of change of pitching-moment coefficient with pitching 
velocity per radian; see appendix E 

rate of change of pitching-moment coefficient with a per 
radian; see appendix E 

pitching-moment coefficient of horizontal tail surface, 
M t/lt S t^t 

acceleration due to gravity, ft/sec/sec 
airplane moment of inertia, slug-ft^ 
airplane radius of gyration about pitching axis, ft 
lift, lb 
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M 

n 


q 

s 

St 

t 

V 

w 




a 

7 

9 

6 

6 


airplane mass, tf/g, slugs 
pitching moment of airplane 
airplane load, factor 

pv2 / 

dynamic pressure, —— , lb/sq ft 

wing area, sq ft 
horizontal-tail area 
time, sec 

true velocity, ft/sec 
airplane weight, lb 

length from center of gravity of airplane to aerodynamic 
center of tail (negative for conventional airplanes), ft 

coefficients of transfer function relating 9 and &; see 
appendix E 

wing angle of attack, radians 
tail angle of attack, radians 
flight -path angle, radians 
angle of pitch, a + 7 
elevator deflection, radians 
downwash angle, radians 



^t 

P 

T 


tail efficiency factor, q^/q 

mass density of air, slugs/cu ft 
dummy variable of integration 
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Matrix notation: 

rectangular matrix 
square matrix 
column matrix 

integrating matrix (see table I) 
rectangular matrices defined in appendix E 

denotes row elements in matrix 
tail 

wing-body combination 

For sign conventions used, see figure 1. 

A dot over a symbol denotes tbe first derivative with respect to 
time, and two dots over a symbol denote tbe second derivative with 
respect to time. 

The symbol A refers to an incremental value. Intermediate 
variables such as An, Al, Acr, Aq>, and Aty and the constant Kpo 
are defined in appendix E. 


I! II 

[] 

{} 
lie 1 1 

IIb|I,||dII,IIeII 

Subscripts: 

i 

t 

WB 


OUTLINE OF METHODS 


The three methods are based on the longitudinal equations of motion 
for horizontal flight and use matrix methods to analyze time histories 
of measured quantities. The equations of motion _used in each of these 
methods are expressed in the form 

+ C L 0® + + C L5 A 5 W 


+ Cm^a + C m g0 + Cmg A5 


— An = C T Ax 
qS 


9 = Cm Ax 
qSc “a. 


( 2 ) 
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These equations apply to a rigid airplane and are based on the 
usual assumptions of linearity, small angles, and no loss in airspeed 
during the maneuver. The equations are further restricted to the range 
in which the variation of the derivatives is linear and also to conven- 
tional wing-tail configurations in which the major contribution to 
damping in pitch is due to the horizontal tail. All the variables are 
given in incremental form measured from a steady- flight trim condition. 

As indicated in reference 2, the four values Ax, a, 0, and A6 
in equations (l) and (2) are linearly dependent; therefore, if four 
simultaneous equations are formed to determine either the force or 
moment derivatives, they cannot be solved uniquely for the unknowns. 

For purposes of analysis the moment equation (2) is integrated 
once and expressed in the form 

■ gZ ® = ^ Ax dt + Ax + A0 + ^ A6 dt (3) 


This form permits the use of numerical integrating methods that are 
more desirable than numerical differentiating schemes when applied to 
flight data. Integrations of the variables are performed by use of an 
integrating matrix ||C|| derived in reference 3 and given in table I 

. ot 

herein. For instance, a time history of / An dt may be obtained 

J 0 

from a time history of An as follows: 



The integrating matrix ||c|| given in table I may be used for any time 
interval At; most of the computations of this paper are based on a 
time inverval of At = 0.1 second. This interval may be too large in 
some cases, and if greater accuracy is desired, a shorter time interval 
may be chosen. 

The essential differences in each of the methods are in the number 
of quantities to be measured. Method A requires four basic measurements 
in time-history form to determine all the derivatives. Method B requires 
three measurements and one supplemental assumption, namely = XC m ^. 

Method C requires two measurements and two supplemental assumptions , 
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x+ 

namely C^ = XC mfo and Cmg = — Cj^. All measurements of flight data 

used are time histories of Incremental values measured from a trimmed 
level- flight initial position. The development of the equations for 
each method is covered in appendixes A to C; in the "body of the report 
the methods are outlined hy stating the pertinent equations in the order 
of computation. Since these computations make extensive use of least- 
squares procedures and are greatly facilitated by the use of matrix 
algebra, most of the equations are given in matrix form. 


Method A 

Of the four basic measurements required with method A, two must 
be incremental elevator angle and incremental tail load and two other 
measurements can be chosen from a possible three, namely incremental 
load factor, pitching velocity, and angle of attack. In this paper, 
incremental load factor and pitching velocity are used. 

The procedure of computation with method A (see appendix A for 
development) is as follows : 

(l) Compute a time history of rate of change of angle of attack a 

from 


a = 0 - | An 


( 5 ) 


or 


{ii} " {h} -f^i} 

P t pt 

(2) Calculate time histories of Aa, A0, / A 9 dt, / AS dt, 

J 0 J 0 

and / / AS dr dt by using the integrating matrix ||c|| and the! 

J 0 J 0 

time histories of a, and AS; for example, 

{aoi} =iic|i^ai^ 


(6) 
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( 3 ) Determine Cj^ and Cpg by least squares from the relation 


Aa + Ci^ A 6 = -^7 An 




qS 


(7) 


or 



-JrH 


(4) Compute the coefficients Kp, K^, K 5 , and Kg of the transfer 

function relating pitching velocity and elevator, deflection by the use 
of the method of reference 4 and the equation 

Kp Ad + K 2 f A0 dt - Kr f AS dt - Kg f f A5 dt dt = -6 ( 8 ) 

Jq J q J n 'J n 


'0 '-'0 


where the measured values of pitching velocity and elevator deflection 
are used. 

(5) Determine Kpo from the relation 




( 9 ) 


by using the results of steps ( 3 ) ana (4). 

( 6 ) Calculate time histories of the intermediate quantities Aqp 
and Ah by using the expressions 


A<p = — 
*t 




( 10 ) 


or 


{"Pi} * sy {*“ 1 } • fa + 0 
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and 


or 


^ = ^ ' ^2 Kl ° ^ ^ Kl0& " Cl S ^ 

■^t 


\j 


(7) Compute Cjh' by least squares from the relation 


( 11 ) 


Cmg A<p = Aj-i 


( 12 ) 


or 


°m, {**} = H 


(8) Determine C^ from the equation 


Cm^ _ ^10 _ 


(13) 


from 


(9) Calculate the time history of the intermediate quantity Act 


Act 


— ~ 9 - Cm* Ax - Cma A0 
qSc ^ u 


m 



dor 
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(10) Compute and Cmg by least squares by using the relation 


C 





Act dt + Cmg 



AS dt = Act 


( 15 ) 


or 


INI 



(ll) Calculate Cj^ and Sj^ from tbe following definitions: 


Cl 6 ~ 5 ^ Cm 0 


(16) 


and 


Cl u ~ ^ Si 


( 17 ) 


(12) Determine the time history of the intermediate quantity Ai|r 
from the equation 


W 


= — An - Cj^B - Cj^a 


( 18 ) 


or 


{**i} = -^{k} 


(13) Compute the refined values of Cj^ and Cl^ by least squares 
by using the relation 


C^ Act + C]Vg AS = At 


( 19 ) 
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or 

These values now include the effects of the Cj,g and Cp^ terms in 
the force equation. 


(llj.) Method A can now he iterated to obtain better values of the 
derivatives by starting the process over at step (5) with the improved 
and Cp^ values from step (13) and following the procedure again. 

The iteration converges rapidly. 



are found from 


( 20 ) 


de . 1 Cm a 

Sa ^ c mg 


( 21 ) 


s - 

35 ti t s t ^ 


( 22 ) 


■' m 0 J 


This procedure shows that 

c c l\ a* 

Cfl \x' C m6 , f ■ J , 


and 


the derivatives Cp^, Cp,£, Cp^, Cp^, Cm^, 
50 Lt 


d5 


may be determined by numerical 


operations on four time histories of measured flight data and through 
the use of the theoretical relationships given as equations (l6), (17)> 
( 20 ), ( 21 ), and (22). 
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Method B 

Three basic measurements are used in method B, one of which must 
be incremental elevator angle and the other two measurements can be 
chosen from a possible three, namely incremental load factor, pitching 
velocity, and angle of attack. In this paper, incremental load factor 
and pitching velocity were used. 

In lieu of the fourth measurement, the supplemental assumption is 
made that 


Cjj^ — XC tTIq 


(23) 


If a value of X is not known in advance, a first approximation (see 
ref. 5) is X = 1. Although the assumption X = ^ imposes a restric- 
tion on the generality of the method, it appears to be justified since 
it reduces computation time to almost one half that required for 
method A and for the examples presented herein gave results which are 
in good agreement with those of method A. 

The method is outlined by merely stating the appropriate equations, 
the development of which is - contained in appendix B. The procedure iB 
as follows: 

(1) Compute the time history of Act. by using equations ( 5 ) and ( 6 ). 

(2) Determine a time history of the intermediate quantity A| from 
the expression 

A| = (1 + X)0 - ^ An (2k) 

or 

{ak} - (1 + M-jV}- - yfaij 



(3) Calculate time histories of 



Ax dt. 


A 6 dt by using the integrating matrix 


of Act, A|, and AS. 



dt, and 


[|C|| and the time histories 
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(10 Compute C^, C m ^, and C m ^ by least squares from the relation 

p t P t p t 

Cn Ax dt + Cm ' Q J A£ dt + Cm§ J AS dt = Q ' (25) 


or 


Sx 

NKCmg^ = ~ <[% 


Cm B 


qSc 


(5) Determine Cj^ from C^ by using equation (23) and then 
determine Cj^ and C]^ by using equations (l6) and (17). 

(6) Calculate the time history of the intermediate quantity Aty 
by inserting these values of C^ and Cj^ into equation (18). 

(7) Compute the values of Cj^ and Cj^ from equation (19). 


(8) The derivatives f • 


/So, 


\da 


i 


dC 


and 


H 


bs 


are then determined from 


equations (20) and (22) and the previously derived quantities. 


Method C 

Method C is an extension of the method presented in reference 4. 
Appendix C contains the development of the pertinent equations upon 
which method C is based. Two basic measurements are used in this 
analysis, one of which must be incremental elevator deflection and the 
other one may be chosen from incremental load factor, pitching velocity, 
angle of attack, and so forth. In this paper incremental pitching 
velocity is used. 

Two supplemental assumptions are made. The first is the relation 
between Cm^ and C m £ given in equation (23) and the second is 


(26) 
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The procedure for method C is as follows: 

(1) Compute the stability coefficients Kp, YLq, and Kg as 

outlined in step (4) of method A. 

(2) Compute Cp^ from the relation 



where 


and 




(1 + X) 



C 2 


(1 + X) 



X K 6 K _ K 6 *t Vl A 

1 + X K^ 1 K 5 1/ 


(3) Determine C^ by using the expression 


(27) 


( 28 ) 


(29) 


I - <l s I Cr 2 + K 1 J 


Cm = - Kp - — _ 
m qSc mV cVm(l + X) 


cVm(l + X) 


'La 


(30) 


(4) Calculate Cmg from 


^ ’ * % ) 


(3D 


( 5 ) Compute C m from 


Cmc — 


I mV 
qSc" qSc" 


Kg 


'La ^ 


(32) 
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(6) Determine Cm^ and Cl§ from equations (23) and (2 6 ) by using 
the values of and Cm§ found in steps (4) and (5), respectively. 

Approximate equations for the stability derivatives are given in 
appendix C. 

In appendix D, method C is modified slightly so that many of the 
stability derivatives can be obtained from a single time history. This 
time history must be the response to an input elevator motion of the 
impulse type. This modified method C comes closest to the ultimate aim 
of this type of analysis, namely to determine the derivatives from a 
single time history. 

One of the important factors in obtaining reliable results with 
the methods outlined herein is the choice of a sufficiently small time 
interval At. In the computations using method C in this paper, in 
one case a time interval of At = 0.1 second was found to be too large 
to give reliable results, and a time interval of At = 0.05 second had 
to be used. 


EXAMPLES 


In order to illustrate the methods outlined in the previous section 
as well as to compare the results obtained, a number of examples are 
given in which the data used are from test runs of a medium jet bomber 
at about the same Mach number. Methods A and C are applied to flights 1 
and 2; whereas all three flights are analyzed by method B. Computations 
are shown in the tables for flight 1, but for the other flights only the 
results are given. 

Table I contains the integrating matrix ||cl| based on Simpson's 
law (ref. 3) which is used in all three methods. 

The airplane characteristics and flight conditions are shown in 
table 11(a) for all three flights. Although the geometric parameters 
are the same, the parameters such as weight, speed, Mach number, center- 
of- gravity position, and altitude vary slightly between the three runs. 

In table 11(b) the coefficients of the transfer function which 
relates pitching velocity to elevator deflection defined by equation (8) 
and computed by the method outlined in reference 4 are shown. These 
preliminary constants are required in methods A and C and the actual 
computations are shown in a subsequent table. 

Time histories of measured and derived quantities for flight 1 are 
shown in table III. The quantities in columns ©, Q), ©, and © are 
measured and the other five quantities are derived from the measured 
quantities. In these tables more decimal places are carried in the 
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measured quantities than are warranted by instrument accuracy in order 
to assure no loss in accuracy in rounding off. The measurements of 
incremental tail load ALt were available only for the times listed, 
and since these covered approximately the natural period of the short- 
period oscillations of the aircraft, the data were considered suffi- 
cient. More of the time histories of the other variables were available 
and were used. 

Method A .- The principal computations illustrating method A are 
presented in table TV; some of the intermediate steps outlined in 
method A are simple computations and are therefore not included in this 
table. Table IV(a) is obtained by applying equation (7) to the data 
given in table III and illustrates step (3) of method A. 

In table IV (b), the computations illustrating the determination 
of Cnig and for steps (7) and (8) of method A are shown. Two 

of the columns are taken from table III and the equations upon which 
the computations are based are (12) and (13). 

Table IV (c) illustrates the computation of and CW for 

step (10) of method A. Two of the columns are obtained by operating 
on columns © and (g) of table III with the integrating matrix ||C|| 
given in table I, and the other column is taken directly from table III. 
The computation is based on equation (15). 

The refined values of and are determined in table IV(d) 

by use of equation (19). Two of the columns are taken directly from 
table III and the other column is derived by use of equation (l8). 

Final results obtained with method A for the data of flights 1 
and 2 after three iterations sire shown in table IV (e). 

Method B . - The principal computations illustrating method B are 
presented in table V. Again, some of the intermediate steps outlined 
in method B are simple computations and are therefore omitted. In 
table V(a) the computation demonstrating the determination of C^, 

C m £, and Cmg by step (!)•) by using the relation (25) is shown. Three 

of the columns are obtained by operating on col umns ©^ © , and © 
of table III with the integrating matrix ||c|| given in table I. 

Table V(b) illustrates step (7), the determination of Ci^ and Cig 

using equation (19)- Two of the columns are obtained directly from 
table III and the other column is derived by using equation (18). 
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In table V(c) final results obtained with method B for three sets 
of flight data are shown. 

Method C .- The principal computations of method C are presented 
in table VI. Table VI (a) shows the computation of K-^, K>), and 

Kg from flight 1 data by the method of reference 4. The integrals in 

table VI (a) were computed by reading the film at 0.05-second intervals 
and using the integrating matrix for At = 0.05 second; this interval 
was necessary in order to obtain reasonable results for the method. 

Use of the time interval At = 0. 1 second did not produce sufficiently 
accurate values of K^ and Kg in this case. Table VI (b) shows the 

computation of the K values for flight 2 data, in this case a time 
interval of At = 0.1 second was sufficiently small to produce reliable 
results for method C. 

In table VI (c) the f ina l results obtained with method C for 
flight 1 and flight 2 data are given along with the results obtained 
by using the approximate formulas of appendix C. 


DISCUSSION 


The three methods presented in this paper are based on the assump- 
tions that the aircraft, has two degrees of freedom (vertical motion and 
pitch), that' the motion of the aircraft can be adequately described by 
the linear differential equations of motion with . constant coefficients 
based on small-perturbation theory, that the aircraft is a rigid body 
with no flexibility, and that the major contribution to the damping 
comes from the horizontal tail. The airplane , its flight condition, 
and the maneuver to be analyzed must therefore fall within the realm 
of these assumptions; that is, the airplane should be operating under 
conditions in the linear range of the coefficients, the maneuver should 
be of the pull-up or push-down variety where little loss in airspeed 
occurs during the maneuver and where displacement angles are small, and 
the maneuver should start from a level-flight trim condition and should 
be in the Mach number range in which these assumptions are valid. 

Since the choice of the method to be used depends primarily on the 
number of measurements which are available, method A is recommended 
when four basic measurements are available, method B when three measure- 
ments are available, and so forth. If, however, an accurate value of \ 
is known in advance, then method B is recommended since it will give the 
most information for the least amount of work. Method C requires more 
work t han method B, and the modified method C is not expected to be so 
reliable as the other methods. 
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In these methods sufficient data to cover the natural period of 
the short-period oscillations of the aircraft should he used. For 
highly damped motions sufficient data should he used to approach the 
steady-state value. 

The accuracy of the results obtained from these methods is influ- 
enced considerably by errors in the instruments and in the record reading. 
Instruments used should be accurate, calibrated both statically and dynam- 
ically, and free from drift and hysteresis. Before an analysis is started 
the data should be corrected for known instrument errors; the records 
should then be read as carefully as possible. Measured tail-load data 
should be corrected for effects of inertia. The accuracy of the analysis 
next depends on the time interval selected for the integrating matrix and 
on the amount of departure from the basic assumptions. Provided the 
initial data are accurate, the smaller the time interval the more accurate 
the results. The differences between the values shown for different 
flights in tables IV(e), V(c), and Vl(c) are believed to represent the 
scatter caused by effects of flexibility, minor nonlinearities, instrument 
errors, record-reading errors, changes in airspeed during the maneuvers, 
and other items which essentially depart from the basic assumptions. 


As may be seen from a comparison of tables IV(a) and IV (d), the 
inclusion of the Cp,^ and Cp^ terms in the force equation for method A 

has little or no effect on Cp^ but has a considerable effect on Cp,g. . 

If the Cp^ and Cp^ terms are retained in the force equation in the 

development of equation (8), the form of. the equation remains the same 
but the K values now include Cp^ and Cp,£ terms. These terms were 


found to have a negligible effect on the K values and their inclusion 
made the equations too unwieldy to handle. For the sake of complete- 
ness, the K values including the Cp^ and Cp^ terms are given in 

appendix E as K x , Kg, i^, K^, K^, and Kg. 


In the actual computation it is recommended that the simultaneous 
equations formed by the least-squares procedure be solved directly by 
the elimination of the variables or by Crout's method. (See ref. 6.) 
The use of a least-squares method permits the calculation of a probable 
error, which is an indication of the fit of the data. The expression 
used in computing the probable error is 



t/5Ii 



MCA TN 2902 


19 


where Bpp is the main diagonal term of the inverted matrix of the 
coefficients , E is the difference between the computed and measured 
value of the variable, N is the number of cases considered in the 
least-squares procedure, and k is the number of unknowns determined 
by the least-squares procedure. , A probable -error analysis was made of 
all the results using flight 1 data and these results are given in 
tables Il(b), IV(e), V(c), and VI (c). This probable-error analysis 
indicates that all the derivatives determined by method A with the 
exception of appear to be more accurate than the derivatives 

determined by method B or C; it also indicates that the derivatives 
determined by methods B and C appear to be of the same order of accuracy. 

When the computed stability derivatives are substituted back into 
the equations of motion, the method that uses the most measurements and 
has the fewest restraints imposed on it would be expected to produce 
the most accurate results and give the best fit to the original data. 

This might not be the case, as illustrated in figure 2 which compares 
the fit of the measured data with the computed data for the three 
methods presented. The results .for method A are more accurate for 
the data herein than the results for method C, but the fit of the 
incremental-pitching-velocity curve for method C is as good as if not 
better than the fit for method A or method B. It appears in general 
that the more coefficients determined from a single time history the 
better will be the fit of the data but the less accurate will be the 
coefficients determined. The fit of the data is interesting since the 
three methods presented are essentially curve-fitting processes in 
which the longitudinal equations of motion are used to fit the flight 
data. A good fit indicates that the equations of motion and assumptions 
used adequately fit the data and the coefficients determined, if inserted 
in the equations of motion, will reproduce the motions of the aircraft. 

In figure 2 the incremental tail load shown for method B was com- 
puted by using the stability derivatives determined from the time 
histories of incremental load factor and pitching velocity. In method C 
the incremental load factor and tail load presented were computed by 
using the derivatives determined from the pitching velocity. These 
time histories indicate how well the derivatives determined on the basis 
of the measurements recorded by one set of instruments will predict thd 
measurements recorded by a different set of instruments. In the case 
of method C the agreement is good; in the case of method B it appears 
that a more realistic value of X than 0.5 should be used. Method B 
is more sensitive to X than method C is. 

Although not presented, the derivatives determined from flight 2 
by methods A and C were used to predict the motions of the aircraft for 
flight 1. A comparison was then made with the actual flight 1 motions 
and it was found that, the predicted motions and the actual motions were 
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in good agreement. These results verify the validity of the method out- 
lined herein as applied to the example airplane. 

A possibility for a further generalized method which would include 
damping effects of wing and fuselage and therefore would make the method 
applicable to the case of swept-wing airplanes may be realized by com- 
bining features of method C with method A in the following manner. Equa- 
tion (A20) may be written in the form 




Now Cmg and may be evaluated by a least-squares procedure, 

t 

provided an accurate value of Cpg is available or can be determined. 
Examination of the probable errors for C^g given in tables IV(e) 
and Vl(c) indicate that, in the case of the medium Jet bomber used in 
the calculations herein, the more accurate value of C;Lg is determined 
by method C by using equations (C12) and (C4) . It is believed this will 
generally be the case for the derivative C^g • It might also be noted 
that this value of CLg will provide more rapid convergence of the 
iterative procedure pf method A. The usual assumption is made that the 
contribution of wing and fuselage to Cjj,^ is negligible. Then 
can be computed through the use of the value of Cm& computed by the 
above procedure and equation (A22) . 

Possibilities for further investigation are to expand the method 
to include flexibility effects and the effects of higher-order deriva- 
tives and to extend the method to the case where the initial conditions 
are known but are not necessarily zero; that is, the maneuvers do not 
start from level-flight trim conditions. The methods could also be 
extended to other configurations such as canard aircraft and tailless 
aircraft, and perhaps a similar analysis could be made of the lateral 
motion of an aircraft to determine the lateral derivatives. 


CONCLUDING REMARKS 


An analysis of longitudinal-stability derivatives by three separate 
methods baa -been presented and applied to flight data. Method A, the 
most general method, requires four measurements in time-history form 
and permits computation of all the longitudinal-stability derivatives; , 
it also requires the most computing time and gives the most accurate 
answers. Method B, which requires three measurements in time-history 
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form and one supplemental assumption, namely = Constant (where 

and C m ^ are the pitching-moment derivatives), gives the most 

information for the least amount of work and gives results which are in 
good agreement with those of method A. Method C requires two measure- 
ments in time-history form and two supplementary assumptions, namely 

— - = Constant and Cmg = — 

Cm 0 c 

effectiveness derivatives, x-fc is the tail length, and c’ is the mean 

aerodynamic chord ). 

The results obtained for the methods presented depend in a large 
measure on accurate instrument measurements and require considerable 
computation to yield adequate engineering answers. Since, however, the 
present trend is toward increased instrument accuracy and expanded 
facilities for machine computation, this direction appears to be the 
one in which flight-data analysis should proceed. 


Cpg (where- Cmg and Cpg are the elevator - 



Langley Aeronautical Laboratory, 
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APPENDIX A 

DETERMINATION OF MNGITUDINAL-STABILITY DERIVATIVES BY USING 
FOUR MEASUREMENTS IN TIME-HISTORY FORM 
Equations of Motion 


The equations of motion for small vertical-plane disturbances may- 
be stated (see fig. 1 for definition) as 





R t qS t Axt + 


&L, 

-§r n t« s t 


(Al) 


[hcS\ _ /SCt\ 5C m , ^Lt- 

19 ' wW 330 + ** + it , ' t9StCt + “§r A8 

(A2) 

where 

Axt = Ax^l - - a 

These equations are for a rigid 

assumptions of linearity, small angles, and no loss in airspeed during 
the maneuver. It should be noted that the variables are all in incre- 
mental form measured from a steady-flight trim condition. 


Substituting equation (A3) into equation (Al) results in 



x t _ q 1 

V &C. V ^ 

body and are based on the usual 
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Since 


An = — 7 
g 


equation (A4) can be expressed as 


(A5) 


where 


W 


— An = Ci^ Aa + Cr^d + C L& 0 + AS 


, T = . (to l\ n St / _ 

^ \da \da /t "k S \ da/ 


(A 6) 


(ATa) 


C I* = 


_(^ l \ _ £t ft 

\da S V 


de 


(A7b) 


CL2 


(^l\ 

\da / 


t " S V ^ 


(A7c) 


Ct _ St 

Cl s B5 n t s 

Substituting equation (A3) into equation (A2) gives 


(A7d) 


) + (&) nt 

St x t 

Axf] 

-h.) 

. x+ de 
- a — S. 

e S. _i_ 

'm / 1 t 

s c 

V 

da/ 

V da 

T 


^t _ 8* ct AC , ^Lfc S t x t 


d5 


TJ + — — A5 + 
Sc 


Tj t — rr AS 
d5 z S c 


(A8) 



2 k 
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or 


10 

qS’c 


= Ax 


lp*\ + Tlt !t _ *£\ 

\da Jy b \da /t S c \ da,/ 


+ a 


T' fj - if 


.fe) 5t af _/ 

\da / t S cv 


: Lt S t x+ 

+ i8| ^r"*Yt + ^r^ b 


,t - (A9) 


which can he expressed as 


10 

— — — Cnw Act + CmA* 1 + CmQ0 + Cmc A5 
qSc 


(AlO) 


where 




= (^Ib + (^ 1"* ^ ^C 1 ■ s) (Alla) 



S t ^ de 

S cV Sa 


(Allb) 



q ?! _i_ 

k B(X / t t S 5y ^ 


(Allc) 


C 


^Lt- 

m 5 = & ^ 


5t £t + £t £t 

Sc Bs lt S c 


(Alia) 


From an examination of equations (All) and (A 7 ) , the following relations 
Eire seen to exist: 



NACA TN 2902 


25 


. . °ma = c Lq (A12b) 

c 

c mg = lfnt ^ Cmg (A12c) 

Equations (A 6 ) and (A10) are linearly dependent in the present 
form and must he suitably altered to be put into a computational form. 
As is well-known, Cl^ and Cp-. sire small, and initially for computa- 
tional purposes the force equation can be expressed as 


l-in = + gjg as 


(A13) 


The derivatives and C-p^ 

tions (A12a) and (A12b) after 


can be determined by means of equa- 
Cjj^ smd Cm^ are determined. 


Equations (A13) and (A10) are now in the identical form of the 
equations of motion developed in reference 7* From figure 1 the 
following relation is seen to exist: 

A9 = Att + A/ (Al4) 

As demonstrated in reference k, equations (A10), (A13), and (Al^) 
may be solved simultaneously to obtain the relation 


A5 + Kg 


which may be expressed in integral form as 


9 + KjQ + K 2 A 9 = K 5 


/ 

J 0 


A6 dt 


(A15) 


9 + K-^ A9 + 


K 2 f A0 dt 



A5 dt + Kg p f A 6 d T dt 

Jo Jo 


where 


Ki = 





(Al 6 a) 
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K 2 ‘ - ^(Sx + C "9 V") 

K 5 = (^6 “ C ^ Cm a) 

Kg = 2§£ qS ^ c _ c c \ 
^ I mV V L a m 6 % “W 


(A l6b) 


(Al6c) 


(Al6d) 


By using "the matrix metbod of reference h 3 Kq^ Kg,, and Kg 

may be evaluated from the time-history measurements of pitching velocity 
and elevator angle. 


Method of Separating and 

This method of separating and Cn^ applies only to conven- 

tional aircraft configurations equipped with a horizontal tail surface 
located to the rear of the wing so that the major contribution to 
damping in pitch is due to the • horizontal tail. 

In order to separate and the tail- load equation is 

developed into a form suitable for computing C m ' Q separately in the 
following manner: 

The incremental tail load is given by 

ALt = ( 3 ^) ntqSt Ax t + nt s tq a 5 ( A 1 7) 


Substituting equation (A3) into equation (A17) gives 


ALt = 


'8C L \ 

^A" t9St 


fi . ii\ . 4 a |i . e h J_ 

\ W V 3a V fa 


AX 1 


6C 


H 


65 


q t Stq AS 


(A18) 
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or 


ALfc 

qS 


/SCt\ 

" V^a j t Tlt 


# - - 



n St BL a 

,Sa y t S V Sa 



S-fc x-t l ^ ^^t S 


T1+ — — -== 9 + 

« c; ir . / 

t 


t 1 ’it 


58 r|t ^ 


which can he expressed as 


ALt 

qS 


cV 


(/^t Cm S _ C m<x) Aa + ^(Cm^a + Cm#®) + c Lg AS 


From equation (Al4) it can he seen that 


Cni^a. + C m^0 — + Cm^ct- + q 7 


hut from equation (Al6a) 


Cm^ + c m§ - 


— Cr - -5= Kn = K. n 
cVm ^ qSc 1 ±u 


Therefore, 


C^a + Cmg0 = Kpo® + Cm£)7 


Substituting equations (A22) and (A23) into equation (A20) gives 


A^t _ 
qS 


" ^2^ Smg 


q - K]_q + C m ^ Ax + — (j^pQa + Cm^T^ + Cpg AS 

u 


(A19) 


(A20) 


(A21) 


(A22) 


(A23) 


(A210 
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which can he expressed as 


f 1 - f? K 10 ^ ^ *10* - 0 ^ AS = 


C . cV 

X * xfc2 


(/^ + 


(A25) 


or 


Afi = Cm^ Acp 


(A26) 


where 


Ap 


= - SL- 


qS x t 


- K 10 Ax - — K 10 a - Cpg AS 


(A27) 


and 


+ *)*<• - + 1)^_ 


(A28) 


From equations (Aik) and (A5) the following relations are self-evident: 


a = 5 - | An 


Ax = A 9 - & [ An dt (A29) 

v J 0 

The relations needed to determine the longitudinal-stability 
derivatives from the flight measurements have now been developed from 
the equations of motion; it remains to express the pertinent relations 
in matrix notation. 


Matrix Form of the Equations 

A powerful tool for data, analysis is provided by matrix methods 
since tabulated time histories .may be conveniently carried in the equa- 
tions of motion. In the matrix solutions using data, it is well-known 
that numerical differentiation is inherently more inaccurate than the 
corresponding integration process. For this reason, whenever necessary, 
the differential equations are expressed in integral form. The first 
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step in matrix solutions is to tabulate the recorded values of the basic 
variables at a number of points to, tp, t 2 , t 3 . . . along a given 
time history as in table III, the interval of time used in most of the 
computations in the paper being At = 0.1 second. These tabulations 
then become the various column matrices ASj_, ALt-p Ani, and 0-^. In 

certain cases smaller time intervals must be used to get reliable results. 
Another means of getting more accuracy is to use integrating matrices 
based on cubic or quartic curves faired through the data in place of 
the parabolic curves. 

The four basic measurements used in the development herein are 
incremental load factor, pitching velocity, tail load, and elevator • 
angle. By use of equation (A29), the time history of incremental 
angle of attack is computed. Equation (A13) may be expressed as 


Ax 0 

A 60 




An 0 

Axp 

ASi 




An-L 

AXg 

• 

• 

A5 2 

• 


l CL 6 

w 

' qS 

An 2 

• 

• 

Aoq 

• 

A 8 n 




• 

Aa n 


<A31) 

Applying least squares, which in matrix notation involves premulti- 
plication of matrix B by its transpose B f , to equation (A3l) yields 



(A32) 
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for which the solution is 



H-b] 


-1 _w_ 

qS 



(A33) 


By the method of reference h, compute K]_, Kqj and Kg from 

the time histories of the pitching velocity and elevator angle. The 
value of K]_q can he obtained from equation (A22). Time histories of 
the derived AM- and Acp functions can now be computed by using equa- 
tions (A27) and (A28) since the value of has been computed from 

equation (A33). 

Equation (A26) becomes 



(A34) 


Applying least squares to equation (A34) results in 




t 

(A^i Acpi) 

0 


(Aq>i ) 2 

0 


(A35) 


From equation (A22) is obtained as 


C n& = ^10 - 


(A3 6) 


Equation (A10) is now expressed in integral form as 


qSc 


9 = C 


nia 


/ 

Jo 


Ax dt + C mA Ax + A 9 + C 


ms 


■'ms 


/ 


AS dt 


(A3 7) 
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which can he rewritten as 


I 

qSc 


® Ax ~ 


Cnig A0 



Now if 


Ac = — j^.z 6 “ Cyn^ Ax — Cmg A0 


then a time history of Acr can he obtained and equation (A 38 ) 
put in the form 


pto 

pto 





/ Ax dt 

J 0 

/ A5 dt 

J 0 




0 

<3 

pti 

pti 





/ . Ax dt 

/ . A6 dt 




Act-l 


J 0 





P t 2 

P t 2 





/ Ax dt . 

J 0 

/ A5 dt 

J 0 


M 

= < 

A°2 

• 

• 

• 

• 


N 


• 

• 

f^n 

/ Ax dt 

J 0 

rt n 

/ A5 dt 

J 0 




• 

* 

Aa n 


»■" fc} • W 

Applying least squares to equation (A4l) gives 


(A3 8) 


(A3 9) 

can he 


(AltO) 


(Altl) 
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and the solution is 



(All-3) 


In order to include the effects of the Cj* and 
initially omitted in the force equation, equation (A 6) 


terms 

is rewritten as 


At = Ax + CLg A5 


where 


Ai}r = An - Ci^9 - Cj^a 


(AMO 


(A^5) 


Method A may now he iterated to obtain more refined values of the 
derivatives. The values of and Cj^ determined by equations (A35 ) > 

(A3 6) , (A 12a) , and (A12b) are inserted into equation (Ah5 ) , and a time 
history of At is computed. New values of and are computed 

from 


(A46) 


or 



If these values of and are used, a new value of K]_q and a 

new time history of Ap can be computed, which, if inserted into equa- 
tions (A34) and (A3 6 ), yield new values of c m£ and C^. The deriva- 
tives Cj^ and CjjQ are again determined from equations (A12a) and 
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(A12b) and a new time history of is computed by using equation (A^5)j 

refined and derivatives are found from equation (A 46). The' 

process converges rapidly. After it has converged, compute and 

Cms from equation (A^3). Thus far C^, C L ^, C^, C^, C^, C^, C^, 

and CniQ have been determined. Then may be determined by 

rewriting equation (Allc) as 

(<Vl\ = _ S5V r • 

^ ) t S t x t 2^ ^ 0 

An examination of equation (A12c) shows that 

de = 5% 

^ Cm 9 

and from equation (A7d) 

c. 

58 ItSt ** 

All the longitudinal-stability derivatives ' are now determined. 1 
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APPENDIX B 

DETERMINATION OF LONGITUDINAL-STABILITY DERIVATIVES BY 
USING THREE MEASUREMENTS IN TIME-HISTORY FORM 
AND ONE SUPPLEMENTAL ASSUMPTION 


The three Basic measurements used in method B are incremental load 
factor, pitching velocity, and elevator angle. The supplemental assump- 
tion made is that Cm^/Cmg is a constant, that is. 


(®l) 

For a first approximation the constant is assumed to he equal to 
1/2 ( see ref. 5) • 

If the definition 


A5 = Xa + 9 = (l + X)0 - X |^An 
is adopted, a time history of A| may he computed. Then 


(B2) 


Cm + Cn^a - Cm^ A| 


(B3) 


The integral form of the moment equation (A37) can then he written 



If time histories of incremental load factor, pitching velocity, 
and elevator angle are measured and ¥, q, S, V, c, and I are known, 
then equation (B4) can he put into matrix form and used directly to 
compute Cju^, C^, and Cmg. 
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The derivative is derived from C m ^ by using equation (Bl) 

and Cq^ and Cq^ are computed from equations (Ai2a) and (A12b) 

C ^u = ^ 

Cl s = Cm e 

These values of Ct. and Cq* are then ins erted into 
equation (A45) ^ 9 

= X. An _ C L - 8 - 

and a time history of At is computed. The values of Cq^ and Cq^ 
are then computed from equation (A44) 

At = Cq^ Ax -t- CLg A5 


Equation (B4) may be expressed in matrix form as 
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or 



Applying least squares yields 



and solving for the derivatives gives 



Equation (A44) expressed in matrix form is 
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Applying least squares to equation (BIO) gives 



and solving for the unkn owns results in 


(BIO) 


(Bll) 



(B12) 


The derivatives 


/SCjA Se 

W )+’ Sot/ 


SC 


and 




"ST 


are determined from equa- 


tions (Allc), (A12c), and (A7d), respectively. 


Thus the derivatives C L a , CLq, CL^, CLg, C m ' Q , C,^ 


(: 


SCl\ Se 


SCi 


Sa > *7 Sa/ 


and 


S6 


are now determined. 
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APPENDIX C 

DETERMINATION OF LONGITUDINAL-STABILITY DERIVATIVES BY USING 
TWO MEASUREMENTS IN TIME -HISTORY FORM AND TWO 
SUPPLEMENTAL ASSUMPTIONS' 


By the matrix method of reference 4, time-history measurements of 
the pitching velocity and incremental elevator angle are used to com- 
pute the K values from the relation 

J. 

8 + K-jQ + K 2 AS = K= AS + Kg r AS dt 

J 0 

(Cl) 

where 


Kl = f ^ f (C* + ' C^) 

(C2a) 

K 2-' T(^ + C "a - ) 

(C2b) 

*5 = 

(C2c) 

*6 " I " C L6 C ma) 

(C2d) 

If it is assumed that 


11 

a* 

0 

(C3) 

and 


Cm 6 = ^ Cl 5 

(C4) 


c 
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these six relations in six unknowns CLg, C^, Cmg, C m ‘ Q) and C^ 

can now be solved simultaneously. The following relations for the 
variables result: 


C Iu= " 



where 


(C5) 



(1 +» 


xtVft 

I 


K X + X ^ 


and 


(c 6) 


c 2 



^ 5 k _ H x tv°A 

1 + XK5 1 K5 I J 


(Cl) 


Cm a = “ 



S§. .... 1 Ct 2 

mV cVm(l + X) 


K X I 


cVm(l + X) ^ 


(C8) 


C m 0 


I 

cVm(l + X) 




(C9) 




°®6 " 


I mV ' Y , 

qSc qSc ^ 


c 


% 

*t 


(CIO) 


Approximate formulas which give a quicker evaluation of the deriva- 
tives with fair accuracy were derived from equations (C5), (C 8 )_, (C 9 ) 
and (CIO) and are 



ho 


NACA TN 2902 




~ mV ^6 
qS K 5 


(Cll) 


Cmg 



(C12) 


C m g 


... l : 1 fe 

1 + X. qSc^5 



(C13) 


(Cll)-) 


The set of approximate formulas has "been found to give results 
which are usually within the accuracy of the method. In table VI (c), 
a comparison is presented between results computed by using the approxi- 
mate relations and the more accurate relations. The set of approximate 
formulas given by equations (Cll) to (Cllf) is used in the development 
of a modified method C which is given in appendix D. 


, ~ I xr. H r I 

~ K 2 - 




k 1 k 6 

(1 + \)K^ 
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APPENDIX D 
MODIFIED METHOD C 


in , 


For some special types of longitudinal maneuvers considerable 
Information may be determined from a single time history. If the 
elevator motion is known to be of the impulse type (a blip of short 
duration) but its magnitude or time history is unknown, then the method 
of appendix C may be modified slightly to yield some of the stability 
derivatives. The method may be used with impulse-type forcing functions 
produced by ballistic devices. If the input is not a pure impulse but 
resembles one, then the modified method may be applied after the 
elevator motion is zero. - Integrals, however, must be evaluated from 
the zero-time trim condition but the least-squares procedure is applied 
only to the time histories after the elevator motion is zero. 

Since the' definite integral of an impulse is a step function and 
the integral of a step function is a ramp function, let • 


f 

J 0 


AS dt = A 


(Dl) 


and 


rr 


AS dT dt = At 


(D2) 


0 u 0 


Substituting these values in the integral form of equation (A15) which 
is 


Pt pt pt PT 

Kq A0 + K 2 / A0 dt - K 5 / AS dt - Kg / / AS dT dt = -0 (D3) 

'-'A 'J Q yj A 'J A 


0 u 0 


results in 


Kq A0 + K 2 f A0 dt - K 5 A - K£At = -0 (D4) 

J 0 
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Equation (iA) may "be expressed in matrix form as 


A 0q 


A01 


A©2 


A0. 


/ 

J 0 


A 0 dt -1 -to 


r* 1 

/ A0 dt -1 -t]_ 


r 


A0 dt 


-1 -t 2 


' n / A0 dt -1 -tn 
J 0 




-01 


Kn 






% 

k^a 




J n 


(D5) 


Equation (D5) is then used to solve for the - stability coefficients K]_, 
Kg, K 5 A, and KgA. 

The following approximate formulas presented in appendix C are 
used to compute the stability derivatives (since only the ratio Kg/Rj 
is used the value of A need not he evaluated) : 


C r » 

^ qS 


(d6) 


°ni0 ~ 


1 I (H _ 
1 + x ^H\ k 5 


K 1 


(D7) 


■' m a 


qSc 


K> I 

!= k 2 - ^ Cm e - — 


% 


qSc 


-Kg - - Jj&f ♦ K ^_ 

1 + X \ K 5/ (! + >0% 


(D8) 
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Also, it is assumed that 




(r»9) 


As indicated previously in appendix A , 



(DIO) 


° L e = c n'e 


(Dll) 


Thus the analysis of a single time history of pitching velocity can 
yield considerable information if it is the response to an elevator 
impulse function; however, the elevator-effectiveness derivatives can- 
not he found by this method. 
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APPENDIX E 

DEFINITIONS OF STABILITY PARAMETERS 


The stability parameters of the methods presented, not previously 
defined in the original list of symbols, can be defined as follows: 




(1 + 




c Lo = 


St Xt 1 

s v.^ 


SCt 

Ci ® * 


St 

s 



St ft 

S C 
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t 


s t x t 2 Be 
S jy Sa 



St x t 2 


^t 

S cV 


an 


Be 


C ®S - 


I*t S-t ^ ^mt 


B 5 


Tit 


S c 


B& 


s t 

Tit g 


K 1 - - t(S + Cm e) 

K 2 = - + c m . %£) 

K5 = ^rls - S Cl sh) 


1(6 " V " Cl 5 C 0 


K 





1+6 
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ALt cV c 

A*i = — p K 10 Ax - — K 10 a - CLc, A5 


qS 


x t‘ 




Ai = (1 + X)6 - X ^ An 
■V 


Act — 9 - Ax - Ciuq A 9 


A<p = — 
x t 


^ 


M - |r Aa - C L jS - O^a 


5, =3|i 


I C Iu 


qSc m. + Ct- 
qS 


mV r . 

„ ^ " ° L e 

" Cm a ^ 


- C 


qS 


+ C 




'. m e 


- "Sx " wfaf 1 U - CL 0 Cm a) 


*2" 


S( u ? ci i) 


qSc 

K 3 - T - 


Cm e c i^ + Cm b(^ " CL e) 


mV 
qS + 
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| C | | is the integrating matrix given in table I. 
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TABLE II.- AIRPLANE CHARACTERISTICS , FLIGHT CONDITIONS 
AND TRANSFER-FUNCTION COEFFICIENTS 
(a) Airplane characteristics and flight conditions 


b, ft 

c, ft 

Center-of-gravity position, 

percent M.A.C 

g/V, l/sec 

l, slug- ft 

k y 2, ft2 

Mach number 

m, slugs 

q, lb/ft 2 -. 

S, ft2 

St, ft 2 

V, ft/sec 

W, lb 

W/qS 

ft 

Ht 

p, slugs/ft3 


Flight 1 

' Flight ' 2 

Flight 3 

89 

89 

8 9 

lit . 016 

14. 016 

14. 016 

27. 3 1 ! 

27.32 

27.44 

0.061923 

0.062854 

0. 062686 

255,865 

255,276 

258,957 

141. 61 

l4l. 61 

141. 61 

0.497 

0.494 

0.496 

1806. 83 

1802.67 

1828. 66 

171 

166 

•171 

1,175 

1,175 

1,175 

289.3 

289-3 

289.3 

520 

512 

514 

58,180 

58,050 

58,880 

0. 289561 

0.297595 

0. 293060 

-33-5 

-33.5 

-33-5 

0.87 

O.87 

0.87 

0.001267 

0.001276 

0. 001281 


(b) Coefficients of airplane transfer function 


Flight 1 


Flight 2 


Coeffic ient 


Probable 


Coefficient 


Probable 


Kl 

4.14 

k 2 

9.547 

k 5 

-9. 767 

K 6 

-14. 624 


It-. 19 

10. 329 
- 10.010 

-15.526 
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.009703 
.055812 
.072886 
. 074625 
. 07 l 4 l 4 
.070698 
.067923 
.068691 
.067923 
.064712 
.043246 
.032077 
.012565 
.022862 
.026023 

.047260 

.075428 

.071623 
. 078674 
.082513 
.088063 
.095707 
.087504 

.030227 

.024956 
.026422 
. 026702 
.019651 
.005864 
.000908 
. 005620 




.03792 

-.12008 

-.35392 

-.62568 

-.92272 

-1.16920 

- 1.40936 

- 1.56736 

-1.74432 
-1.85808 
- 1. 75064 
- 1.75064 
- 1.45360 
- 1.25768 
-.88480 
-.46136 
-.05056 
.28440 
. 72680 
1.17552 
1.46624 
1.73168 
2.07296 
1. 85808 

1.78224 
1. 62424 
1.46624 
1.22608 
.99856 
.75208 
.34760 
.29072 
.01264 
-.18960 


0 

0 

-.019808 

-.074064 

-.125305 
-. 149419 
-.165781 

-. 167504 
-.163628 
-.155447 
-. x 42960 
-.127027 
-.090426 
-.053810 
.014640 
.058562 
.091287 
.133055 
.167503 
.196353 
.207549 
. 214008 
.217022 
.217022 
.186880 
.121859 
.051241 
.006889 

-.039516 

-.072772 

-.080953 

-.088704 

-.083537 

-.085259 

-.068896 

-.062437 


-.000922 
-.005175 
-.013921 
-.024871 
-.035932 
-.046216 
— 054782 
-.061513 
-.066176 
-.068518 
-.068301 
-.064585 
-.056328 
-.044150 

-.030073 
-.014737 
.001884 
. 019360 
. 036467 ! 

.051697 

.065082 

.076895 

. 085314 

.088583 

•.085705 

.077804 

.066452 

.052350 

.037762 

.023846 

.012066 

.001887 

-.006638 

-.012536 


0 

0 

-.030886 
-. 107378 
-.177000 
-. 204757 
-.220103 
-. 215056 
-.201806 
-.184643 
-.160433 
133012 
-.081436 
-.026512 
. 066966 
. 126783 

. 164326 
.213867 
. 252820 
.285725 
.288821 
.284616 
. 280136 
. 271918 
. 216138 
.125260 
. 021681 
-. 039955 
-. 104671 
-. 147119 
-. 152347 
-. 156342 
-.136064 
-. 136890 
-.103735 
-. 087786 


.010591 

.068067 

. 183905 
.328389 
.474922 

. 610391 

. 724134 
. 812714 . 

.875829 

.908173 

.902665 

.856022 

.744686 

.586752 

.400397 

.196931 

-.022338 

-.250374 

-.476562 

-.679353 

-.854892 

- 1.010045 

- 1.124562 


0 

. 000741 
. 004929 
-.003122 
-.022667 
-.046165 
-.070393 
-.092683 
-.112697 

-.127257 
-.138772 
-.148055 
-. 147554 
-.145166 
-.125729 
-. 100705 

-.072582 

-.039441 

-.005166 

.031245 

.067093 

. 100087 
.127688 
. 156658 
.183397 



.025122 

.028773 

.031746 

.034084 

.035700 

.036590 

.034981 

.032478 

.025852 

.020274 
.014992 
.008194 
. 00133 ^ 
.005589 
.011195 
.016361 
. 021140 
.025542 
.026734 
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Ax 

(table in. 
column (6) ) 

A5 

(table HI. 
column (2) ) 

^-An 

qS 

0 

0 

0 

0 

.009703 

0 

-.000922 

.055812 

. OIO98O 

-.005175 

.072880 

-. 034770 

-.013921 

. 074625 

-.102481 

-.024871 

. 071414 

-.I8n72 

-.035932 

. 070698 

-.267184 

-.046216 

.067923 

-.338554 

-.054782 

. 068691 

-.408095 

-.061513 

. 067923 

-.453846 

-. 066176 

. 064712 

-.505087 

-.068518 

.043246 

-.538027 

-.068301 

. 032077 

-.506917 

-.064585 

-. 012565 

-.506917 

-.056328 

-.022862 

-. 420906 

-.044150 

-. 026023 

-.364175 

-.030073 

-.047260 

-.256203 

-.014737 

-.075428 

-.133592 


Ax 

(table III, 
column (6) ) 

' 

A5 

(table III. 
column (§) ) 

. ^An 

qs 

0.001884 

-0. 071623 

-o.oi464o 

. 019360 

-.078674 

.082351 

. 036467 

-.082513 

.210453 

. 051697 

-.088063 

.340384 

. 065082 

-.095707 

.424566 

. 076895 

-.087504 

.501427 

. 085314 

-. 030227 

.600248 

. 088583 

-.024956 

.538027 

. 085705 

-.026422 

.516067 

.077804 

-.026702 

.470316 

. 066452 

-.019651 

.424566 

.052350 

-.005864 

.355025 

.037762 

.000908 

.289144 

.023846 

. 005620 

.217773 

. 012066 

.007958 

. 100651 

.001887 

. 010611 

. 084181 

-.006638 

. 009i45 

.003660 

-.012536 

. 010175 

-.054901 


CI^ Ax + CLg AB = — An 


0.087761 -0. 0533^5] (CiJ 

0.053345 0.100964_ [CLg 


°Iu = 7. °T 


0]^ = 0.262 


\NACA, 
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TABLE IV. -• COMPUTATIONS ILLUSTRATING METHOD A - Continued 
(b) Determination of C^ ami C^ by steps (7) and (8) 


t 

Acp 

(table III, 
column (5) ) 

Ap 

(table III, 
column (^) ) 

t 

Ap 

(table III, 
column (8) ) 

Ap 

(table III. 
column @ ) 

0 

0 

0 

1.2 

0. 902665 

-0. 147554 

.1 

0 

. 000741 

1.3 

. 856022 

-.145166 

.2 

. 010591 

. 004929 

1.4 

.744686 

-.125729 

.3 

. 068067 

-.003122 

1.5 

. 586752 

-.100705 

.4 

. 183905 ’ 

- . 022667 

1.6 

. 400397 

-.072582 

.5 

. 328389 

-.046165 

1.7 

- 196931 

-.039441 

.6 

.474922 

-.070393 

1.8 

-.022338 

-.005166 

.7 

.610391 

-. 092683 

1.9 

-.250374 

. 031245 

.8 

. 724134 

-.112697 

2.0 

-.476562 

. 067093 

• 9 

. 812714 

-. 127257 

2.1 

-.679353 

. 100087 

1.0 

. 875829 

-. 138772 

2.2 

-.854892 

. 127688 

IBS 

. 908173 

-.148055 

2.3 

-1.010045 

. 156658 

m 



2.4 

-1. 124562 

. 183397 


Cmg |Acpi| = 


9-9342640^ = -1.578836 


Cmg = -0.1589 


Cmi - K]_q - Cma 


Cm^ = -0.0803 
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r Ax dt 
J 0 

( A5 dt 
J 0 

Act 

0 

0 . 

0 

0 

. 000182 

0 

-.000018 

'. 003154 

-.001986 

-.000295 

. 009716 

-.007957 

-.001232 

.017219 

-.014935 

-.003153 

. 024501 

-.020224 

-.006200 

.031585 

-.025122 

-. 010314 

.038487 

-.028773 

-.015379 

.045288 

-.031746 

-.021209 

.052139 

-.034084 

-.027613 

. 058791 

-.035700 

-.034367 

. 064103 

-.036590 



f Ax dt 
J 0 

[* A8 dt 
J 0 

Act 

-0.041237 

0. 067784 

-0.034981 

-.047910 

. 068473 

-.032478 

-.053989 

. 066415 

-.025852 

-.059045 

. 064119 

-.020274 

-.062767 

. 060601 

-.014992 

-. 065018 

. 054201 

-.008194 

-.065668 

. 046582 

-.001334 

-.064612 

. 039040 

. OO5589 

-.061805 

. 030954 

.011195 

-.057382 

.022443 

.016361 

-.051530 

. 013272 

. 021140 

-.044418 

.003702 

.025542 

-.036264 

-.002593 

.026734 


Ax dt + C n 


A& dt 


0. 042893 -0. 036635 JCn^ 

0.036635 0. Oi)-5583 [Cmg 


0.005959] 

-0.018185 


= -0. 644 


Cmg = -O.916 


















56 


NACA TN 2902 




Aa 

(table III. 
column © ) 

A5 

(table HI. 
column (2 ) ) 

A\Jr 

0 

0 

0 

0 

.009703 

0 

-.000922 

.055812 

.013041 

-.005175 

.072880 

-.027606 

-.013921 

.074625 

-.090675 

-.024871 

.071414 

-. 167519 

-.035932 

. 070698 

-.'252511 

-.046216 

.067923 

-. 324221 

-.054782 

. 068691 

-.394650 

-.061513 

. 067923 

-.44154-9 

-.066176 

.064712 

-.494407 

-.068518 

.043246 

-.529178 

-.068301 

.032077 

-.501508 

-.064585 

-.012565 

-.505173 

-.056328 

-.022862 

-.425395 

-.044150 

-.026023 

-372653 

-.030073 

-.047260 

-.267181 

-.014737 

-. 075428 

-. 147869 



Aa 

(table HI. 
column © ) 

AS 

(table III. 
column (2) ) 

A+ 

0.001884 

-0. 071623 

-O.O315H 

. 019360 

. -.078674 

. 063290 

. 036467 

-.082513 

.191190 

.051697 

-. 088063 

. 321409 

. 065082 

-.095707 

.405894 

. 076895 

-. 087504 

.483307 

. 085314 

-.030227 

.585856 

. 088583 

-.024956 

.529696 

. 085705 

-.026422 

.514646 

.077804 

-.026702 

.473006 

.066452 

-.019651 

.431572 

. 052350 

-.005864 

.364860 

.037762 

. 000908 

. 299324 

.023846 

. OO562O 

.228216 

.012066 

.007958 

.109736 

.001887 

. 0106ll 

.093319 

-.006638 

.009145 

. 010582 

-.012536 

. 010175 

-.049055 


Aty = Gi^ Aa + Gig A6 


0.087761 -0.053345 JCLo, 

0. 053345 0. 100964 IcLg 


0.597599 

-0.331104 


C I U - T-09 


Ct. . 0.469 
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TABLE IV.- COMPUTATIONS ILLUSTRATING METHOD A - Concluded 
(e) Final results using method A after three, iterations 


Cm 6 

(dCi/da) t 
Se/Sa 
■ SCL 


■ Flight 1 

Probable 
error for 
flight 1 

Flight 2 

7.09 

, °- 1]L 3 

6.70 

0.468 

0. 105 . 

0.446 

0.072 

0.0004 

0.076 

0.028 

0.0004 

0.033 

-0. 622 

0.003 

- 0 . 711 

-0.171 

0.001 

- 0 . 181 

-0.068 

0.001 

-0. 078 

-0.914 

0.003 

-O.968 

4.84 

> 

5-05 

0.425 


o:46i 

2.19 


2.08 
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TABLE V. - CCMPDTATIOHS ILLUSTRATIHG METHOD B 
(a) Determination of Cn^, Cmg, and 


1 

/ Ax dt 

f Ag dt 
J 0 

/ A5 dt 

■ 

1 

Ax dt 

m 

/ A5 dt 

-S-& 

1 

4 ) 

J 0 


H 

J 0 

Hi 

J 0 

qSc 

0 

0 

0 

0 

0 

1.8 

-O.065668 

-0. 104508 

0.046582 

0.015219 

.1 

0 

0 

.000182 

0 

1.9 

-.064612 

-.077530 

.03901*0 

.017840 

.2 

-.000018 

-.001164 

.003154 

-.001800 

2.0 

-.061805 

-.048742 

.030954 

.018857 

.3 

-.000295 

-.007697 

.009716 

-.006729 

2.1 

-.057382 

-.020009 

.022443 

.019444 

.4 

-.001232 

-.022265 

.017219 

-.011385 

2.2 

-.051530 

.008260 

.013272 

.019718 

-5 

-.003153 

-.041702 

. 024501 

-.013576 

2.3 

-.044418 

.035893 

.003702 

.019718 

.6 

-.006200 

-.063115 

.031585 

-.015062 

2.4 

-.036264 

.060589 

-.002593 

. 016979 

• 7 

-.010314 

-.085043 

.038487 

-.015219 

2.5 

-.027526 

.077951 

-.005296 

.011072 

.8 

-.015379 

-.105918 

.045288 

-.014867 

2.6 

-.018770 

.084949 

-. 007809 

.004656 

-9 

-.021209 

-.125274 

.052139 

-.014123 

2.7 

-.010553 

.083685 

-.010526 

.000626 

1.0 

-.027613 

-.142554 

.058791 

-.012989 

2.8 

-.003317 

. 076269 

-.012905 

-.003590 

1.1 

-.034367 

-157253 

.064103 

-.011541 

2.9 

.00261*6 

.063493 

-.014122 

-.006612 

1.2 

-.041237 

-.168003 

.067784 

-.000216 

3.0 

.007146 

.048510 

-.014312 

-.007355 

1.3 

-.047910 

-.173429 

.068473 

-.004889 

3.1 

.010221 

.033065 

-.013966 

-.008059 

1.4 

-053989 

-.171125 

.066415 

.001330 

3.2 

.012003 

.018620 

-.013267 

-.007590 

1-5 

-.059045 

-.161157 

.064119 

.005321 

3.3 

.012687 

. 005148 

-.012304 

-.007746 

1.6 

-.062767 

-.146702 

.060601 

.008294 

3.4 

.0121*28 

-.006739 

-.011282 

-.000280 

1.7 

-.065018 

-.127892 

.054201 

.012089 

3.5 

.011447 

-.016172 

-.010294 

-.000869 




pt 


ot 


nt 





- 4:6 = 

qSc 

c % / 
oq 

dt + Cttiq 

X 

dt + Cnie / AS dt 

J o 



0.0448786 0.0743828 -0.0370442 

0.0743828 0.2959167 -0. 1121529 

-0.0370442 -0. 1121529 0.0471095 


Cm^ — -0. 624 



.00720991 
0.0059555 
1-0.0006984 I 


Cmg = -0.149 


Cmg = -0.861 
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TABLE V.- COMPUTATIONS ILLUSTRATING METHOD B - Continued 
(b) Determination of Cj^ and CLg 


B 

Ax 

(table III. 
column (§) ) 

A5 

(table III, 
column (2) ) 

At 

1 

Ax 

(table III. 
column (§) ) 

AS 

(table III. 
column (2) ) 

Aijr 

0 

0 

0 

0 

1.8 

0.001884 

-0.071623 

-0. 030434 

.1 

0 

.009703 

0 

1.9 

. 019360 

-. 078674 

.064502 

.2 

-.000922 

.055812 

.012909 

2.0 

. 036467 

-.082513 

. 192411 

.3 


. 072880 

-.028063 

2.1 

.051697 

-. 088063 

.322604 

.4 

-.013921 

. 074625 

-. 091424 

2.2 

. 065082 

-.095707 

.407066 

.5 

-.024871 

.071414 

-.168381 

2.3 

. 076895 

-.087504 

.484441 

.6 

-.035932 

. 070698 

-.253435 

2.4 

. 085314 

-.030227 

.586746 

.7 

-.046216 

. 067923 

-.325120 

2.5 

. 088583 

-.024956 

. 530203 

.8 

-.054782 

.068691 

-.395488 

2.6 

.085705 

-.026422 

.514713 

• 9 

-.061513 

.067923 

-.442312 

2.7 

.077804 

-. 026702 

.472812 

1.0 

-.066176 

.064712 

-.495065 

2.8 

.066452 

-.019651 

.431104 

l.l 

-.068518 

.043246 

-.529718 

2.9 

. 052350 

-.005864 

. 364215 

1.2 

-.068301 

.032077 

-.501830 

3.0 

.037762 

.000908 

.298661 

1.3 

-.064585 

-.012565 

-.505261 

3.1 

. 023846 

.005620 

.227539 

1.4 

-.056328 

-.022862 

-.425090 

3.2 

. 012066 

.007958 

. 109150 

1.5 

-.044150 

-. 026023 

-.372095 

3.3 

.001887 

. 010611 

.092732 

1.6 

-.030073 

-. 047260 

-.266469 

3.4 

-.006638 

.009145 

.010140 

1.7 

-.014737 

-.075428 

-. 146952 

3.5 

-.012536 

. 010175 

-.049418 


Ai}r = Aa + CL5 AS 


0. 087761 
-0.053345 


Clu = 7.0 9 


-0.053345 

0.100964 



0.598132| 
1-0. 332334 | 


Clq = O.456 
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TABLE V.- COMPUTATIONS ILLUSTRATING METHOD B - Concluded 

(c) Final results from three sets of flight data using 
method B with X = 0.5 



Flight 1 

Probable 
error for 
flight 1 

Flight 2 

Flight 3 

C I« 

7.09 

0.113 

6.68 

6.78 

% 

0.1*56 

0.105 

0.1*02 

0.1*20 


0.062 

. 0.008 

- 0.057 

0.053 

°Id 

0. 031 

O.OOl* 

0.028 

0. 026 

Cm a 

-0. 62h 

0.026 

-0. 670 

-0. 698 

Cjn 0 

-0. ll*9 

0.019 

-0. 136 

-0. 126 

Cma, 

-0.075 

0.010 

-0.068 

-0.063 

Cm 6 

-0.861 

O.O63 

-0. 813 

-0.892 

(SCL/Sa) t 

1*. 21 


3.79 

3.52 

^CLt/^5 

2.13 


1.87 

1.96 






(a) r©t«rminatloo of X\j *£> *6 flight 1 data 


pt 


-Jo 


At dt 


pt pT 

-/ / Afl dT dt 

\JQ '-'0 


-A0 




pt 

pt 

pt pT 

t 

AS 

/ AS dt 

-/ A5 dt 

-/ / Aft dr 



Jo 


J 0 

1.2 

i-3 


-0.062600 

-.076111 

-0.06933 

-.07124 

-.06946 

-0.03663 

-.04570 

1.1 

- , 142^56 

-.090671 

-.03271 

1.5 

-.130760 

-.101769 

-.06726 

-.05958 

1 £ 

=. 1*1333 

_ 1 1 ftOfM 


_ oC/iifi 

u 

-.120101 

*.104535 

-.1W0T 
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TAM n.- CCKPtHAnore HimTUATlX HEM 0 - CoctintMd 
(b) DotsrsJLnatloc of K^, Kg/ Kj, ud E^ from flight 2 data 
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TABLE VI.- COMPUTATIONS ILLUSTRATING METHOD C - Concluded 
(c) Final results using method C with X = 0.5 
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Figure 1.- Sign conventions employed. 
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